On Shear-Free perturbations of FLRW Universes 
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A surprising exact result for the Einstein Field Equations is that if pressure-free matter is moving 
in a shear-free way, then it must be either expansion-free or rotation-free. It has been suggested this 
result is also true for any barotropic perfect fluid, but a proof has remained elusive. We consider 
the case of barotropic perfect fluid solutions linearized about a Robertson- Walker geometry, and 
prove that the result remains true except for the case of a specific highly non-linear equation of 
state. We argue that this equation of state is non-physical, and hence the result is true in the 
linearized case for all physically realistic barotropic perfect fluids. This result, which is not true 
in Newtonian cosmology, demonstrates that the linearized solutions, believed to result in standard 
local Newtonian theory, do not always give the usual behaviour of Newtonian solutions. 



I. INTRODUCTION 

This paper deals with a number of interesting proper- 
ties of shear-free perfect fluid solutions of General Rel- 
ativity (GR). The motivation for this work stems from 
the desire to probe the relationship between relativis- 
tic and Newtonian cosmology and their implications for 
the study of the growth of large-scale structure in the 
Universe. Of particular importance is understanding the 
differential properties of time-like geodesies which de- 
scribe the fluid flow in cosmology. The kinematics of 
such fluid flows are described by the expansion G, shear 
(or distortion) aab, rotation w'^, and acceleration Aa of 
the four-velocity field u" tangent to the fluid flow lines. 
Their governing equations are obtained by contracting 
the Ricci identities (applied to m°) along and orthogonal 
to u", which determine how they couple to gravity via 
the Einstein Field Equations [l| . 

Of particular interest is what role the shear plays in the 
relationship between Newtonian and relativistic cosmolo- 
gies. For example it has been known for some time that 
quasi- Newtonian descriptions of cosmology, the so-called 
Silent models, may be constructed for observers which 
move along geodesies which are both shear-free and irro- 
tational 2] . The intricate relationship between the kine- 
matic quantities in Newtonian and relativistic cosmolo- 
gies is most strikingly seen in a remarkable result first ob- 
tained by one of us in 1967 Q. In this paper it was found 
that if the four velocity vector field of a barotropic perfect 
fluid with vanishing pressure is shear free, then either the 
expansion or the rotation of the fluid vanishes. This is a 
purely local result to which no corresponding Newtonian 
equivalent appears to hold, as counter-examples can be 
explicitly constructed [J]. Given that this theorem and 
its extensions appear to hold for arbitrarily weak fields 
and for fluids of arbitrarily low density, one needs to un- 
derstand why the Newtonian approximation fails. 



The result has be extended to general barotropic fluids 
for number of special cases by Senovilla [5j, but has yet 
to be proved in general. As a first step towards this goal, 
we examine what ever result holds in situations where 
the hydrodynamic and gravitational equations have been 
linearised. Of course there are many ways of doing this, 
but one way that is cosmologically relevant is to linearise 
the equations about a Friedmann-Lemaitre-Robertson- 
Walker (FLRW) background ^^. These almost FLRW 
models can be thought of as lying somewhere between the 
full non-linear GR situation and Newtonian theory, at 
least in the cosmological context, and therefore an anal- 
ysis of theorem in this context could shed some light on 
the generality of the result. We show that it remains true 
for such linearised barotropic perfect fluid solutions, un- 
less the fluid obeys a highly non-linear equation of state 
(see (j55p below) which we argue is non-physical. Hence 
the result remains true for physically realistic equations 
of state in an almost-FLRW geometry. 

This result will be useful in obtaining and studying 
new perfect-fluid solutions of Einstein's field equations 
with a shear-free velocity vector field, and in examining 
how linearized General Relativity solutions relate to the 
Newtonian case, which is the foundation of astrophysical 
studies in cosmology. 



II. LINEARISED FIELD EQUATIONS ABOUT 
FLRW BACKGROUND 



To perturb the FLRW spacetime we use the standard 
1-1-3 covariant approach ilj], where we must first define a 
time- like congruence with a unit tangent vector m°. The 
natural choice of this vector is tangent to the the matter 
flow lines. Then the spacetime is split locally in the form 
R ^ V where R denotes the worldline along u°- and V 
is the 3-space perpendicular to u"-. Then any vector X'^ 
can be projected on the 3-space by the projection tensor 
h°'t, = g°'t, + u°'ui,. At this point, two derivatives are 

deflned: the vector u"" is used to deflne the covariant 
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by a dot) for any tensor T°" C.d, given by 

rpa..h ^,G-\1 rpa..h 

^ c.d — U Vgi c.d 



(1) 



and the tensor hat is used to define the fully orthogonally 
projected covariant derivative D for any tensor T'^'^c.d- 
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with total projection on all the free indices. Angle brack- 
ets denote orthogonal projections of vectors, and the or- 
thogonally projected symmetric trace-free PSTF part of 
tensors: 
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This splitting of spacetime also naturally defines the 3- 
volume element 



eabc = -V\9K5l5l5%u\ 
with the following identities 
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The covariant derivative of the time-like vector u"" can 
now be decomposed into the irreducible parts as 

Vo-Ub = -AaUi, + -habQ + (Jab + <^abc^'' , (6) 

where Aa = via is the acceleration, 9 — Dau'^ is the 
expansion, aab = Di^a'^b) is the shear tensor and w"" = 
e'^'"^Di,Uc is the vorticity vector. Similarly the Weyl cur- 
vature tensor can be decomposed irreducibly into the 
Gravito-Electric and Gravito-Magnetic parts as 



Eab = CabcdU^U 



E. 
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(7) 
which allows a covariant description of tidal forces and 
gravitational radiation. 

In the 1-1-3 covariant perturbation theory [6l-[l3||. we 
consider the background to be FLRW where the Hubble 
scale sets the scale for the perturbations. The quantities 
that vanish in the background spacetime are considered 
to be first order and are automatically gau ge-i nvariant 
by virtue of the Stewart and Walker lemma [ij| . In the 
perturbed spacetime the matter is considered to be a 
perfect fluid with the Energy Momentum tensor 



T"" = {fi+p)u''u''+pg 



^ab 
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so the vector field u°' is uniquely defined as the timelike 
eigenvector of T°-^ as long as ^ -I- p 7^ (the heat flux 
g" and the anisotropic stress TTab vanish in the perturbed 
spacetime). Furthermore, we assume the matter to have 
a barotropic equation of state p = pifJ-) satisfying the 
Weak and Dominant energy conditions. We exclude the 
vacuum case, therefore the energy conditions will be 



for both the background spacetime and the perturbed 
solution (and the Minkowski and De Sitter backgrounds 
will not occur). The local isentropic sound speed is 



C2 = ^ . < C^ < 1 



(10) 



The bound on the local sound speed is required for local 
stability of matter (lower bound) and causality (upper 
bound), respectively. 

Now we consider shear-free perturbations and hence 
the shear tensor {(Jab) vanishes identically. With the con- 
ditions above, the linearised field equations are then as 
follows: 



Propagation equations 



^l>0 ; ^l+p>0 ; M> IpI 
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Constraint equations 



(Co)"^ := E"'' ~ D^^A^^ = , 



[CiY := £>°e - -e^^^DbUJc = 



(C2) := D-LJa = , 



(Ca)"'' -.^ H"'' + D''''lo''^ =0 



{CUr:^D-p+{p+p)A''=0, 



{C5T — DbE"'' - -D^fi = 



{CeT -.^ DbH'^' + Oi + p)uj'' ^ 
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(17) 
(18) 
(19) 
(20) 

(21) 
(22) 



We note that the constraints (Ci)", (C2), (C3)"^ (C^T 
and {Cq)"" are the constraints of the Einstein field equa- 
tions for general matter motion specialized to the shear- 
free case and are known to be consistently time propa- 
gated along M° locally. However the conditions aab — 
and (7° = give the two new constraints (Cq)"^ and (C4)'' 
respectively. 

We also use the following linearised commutation rela- 
tions for shear-free congruences: For any scalar '/' 



Using equation (l30l) in (|29|) we get the constraint 



[DaDb - DbDa]f = 2€atcU}^f , 



(23) 



If the gradient of the scalar is of the first order, wc then 
have 

[D'^DbDa - DbD^]f = \U- \q^^ Dbf , (24) 



+2edbcD''iLu'f) . (25) 

Also for any first order 3- vector V^ = 1^^"^ , we have 

2^ l"2\,a. 



[D-Db - DbD-]Va =3^-36') h'^iaVb] , (26) 



h\h\{DdV') = DbV(-) - -QDbV"^ (27) 



h" ciD'^V) ^ Db{D^^V°-^) - -QD'^V . (28) 

o 

Using the field equations and identities of this section 
we will now investigate the compatibility of the new con- 
straints with the existing ones in terms of the consistency 
up to the linear order of their spatial and temporal prop- 
agation. 
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For the new constraints (Co)"'' and (C4)" to be compat- 
ible with the existing ones, the constraint {Ci)b must be 
satisfied. 

To check the spatial consistency of {C'j)b on any initial 
hypersurface we take the curl of ([51]) to get 



r^ 1 

s ^acb i-i Fjc^ Fi 7~i ,, I ^acb 1 



e'^^'D^D^'D^^Da))! + -e'^'^^D^Dbpi - , (32) 
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which using (P5)) gives 

^^acb^^^rf^ ^^ + ^c^> = . (33) 

^Ji + p 3 

Breaking the PSTF part according to equation ([3]) and 
using the commutators ([M|) . (^51) we have 



-^s neb 



fl + p 



'^D.DbD^fi +1^- l^A DcDbpi 

+ fiedbkOcD'^Lu''] + ^cj"/i = .(34) 



Again using ((23|) and (O in the above equation we get 

H + p [3 V 37 

+ fiD^Lj''] + |w"/i = . (35) 

Now from the relation ([25]) and using ^TE\\ we know 

DkD^oj'^ = l{^^- J©') ^" , (36) 

Plugging ([55]) and p3|) in ([55]) and simplifying wc finally 
get 



III. CONSISTENCY OF THE NEW 
CONSTRAINTS 

The conditions of shear-free perturbations and the 
matter being a perfect fluid in the perturbed spacetime 
give the new constraints (Co)"** and (C4)" respectively. 
To check their compatibility with the linearised existing 
constraints of Einstein field equations (henceforth all the 
equations are up to the linear order), we plug (Co)bd in 
(6*5)6 to get 



D'^Di^bAd) - \Db^i - . 

Now from the constraint {Ci)b we have 
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{CsT := 



-u^Y + ciD^Lj" 



0, 



where 



Y = ^, + p + ci[^l--Q^ 
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(38) 



From (Cg)" we can immediately see that for mat- 



ter with constant pressure (p = constant 



0), 



shear-free perturbations are consistent iff 80;" = (as 
according to the second condition of ([9]), /i -f p > 0). 
That is, if the geodesies of the matter congruence in 
the perturbed spacetime are shear-free then they should 
be either expansion- free or vorticity-free (or both). This 
shows that the results of \^] and |J/ for pressure- free mat- 
ter are true for the linearized theory. 

However for a general equation of state, all we can say 
from the equation (|37p is, either the matter congruence 



is expansion free {Q = 0), or the vorticity vector must 
satisfy 






(39) 



for the new constraints to be spatially consistent on any 
initial hypersurface. 

Now let us check the temporal consistency of the con- 
straint (p9|. Propagating it along u'^ we get 



We can easily see that 

c2 = -e(^ + p)0. 

Now from (|28p we have 

We know from the constraint (|18p that 
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Using (P7)) the equation (^5]) becomes 



DbiD'^'cj-^)- ^ -Dt 



D^uj(^) - -QD^'uo'' 



+D"ujlb) - -QD^ijj^ 



(44) 



Simplifying the above equation using p2|) , (|20)) and (l23t , 
we get 

DbiD^^'u;''^} = -^0(1 - c^)(D2a;" + DbD^u'') . (45) 

Putting equation psj) in (jH]), we have 

where 



c? 5 



2 6 
Using (123) and dMl, (SSI) becomes 
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Combining (HTI) and pS)) and using ([5^ we get 
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Also from dTT]), dTH]) and (gT]) we have 
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(m+p)(1 + c2) + 



2 2, 

3C. I M 



^«^ 



(51) 



Now using dn]), ([201), ([231) and (l50|) we get 



l{u;-Y) = -'^-u'^Q 
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Finally using (|49|) and ([52|) in (|40| and simplifying, we 
get 
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where ^R = 2[/i — (1/3)0^] is the spatial curvature. In 
FLRW spacetimes it can be written in term of the scale 
factor 'a(i)' as, 
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where k = —1,0,-1-1 denotes open, flat and closed uni- 
verses respectively. Thus we can easily see that for the 
new constraints to be spatially and temporaly consistent 
we must have either a;"0 = or the barotropic equation 
of state must satisfy the following non-linear higher order 
DE: 
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(55) 



We see that the shear-free results of |J/ and fj/ are 
avoided, at least at the linearised level, if the equation of 
state of the matter solves if?3|). However, a priori it seems 
highly unlikely that any realistic barotropic equation of 
state will obey this extremely non-linear equation. We 
now try to find solutions of this equation, under various 
simplified assumptions or realistic initial conditions, to 
confirm it is nonphysical. 

1. Flat universe (fc — 0) with c^ = constant 7^ 0: 

This is the simplest case in which the equation ((55)) 
reduces to a simple algebraic equation 



5 1 . 

6 + r^ 



=0 



(56) 



which gives c^ = —5/3. This is physically not pos- 
sible as the lower bound on the local sound speed 
(fTUI) is violated. 



2. Closed/open universe with c^ = constant ^ 0: In 



this case also, the equation ([55|) reduces to an 
algebraic equation, and we get the relation 



(57) 






Differentiating ([57)1 with respect to /i and using ([2 
we get 



2 ^i? ^ _2.(| + i^^) ., , „2^ 



-(1 + c^) 



(58) 



3(/z+p) c2(l-c2) 

Eliminating ^R/{^ + p) from (|57l) and ([55)) we get 
the solution c^ = —1/3, which again violates the 
lower bound of the local sound speed. 

3. Flat universe with varying sound speed: In this 
case the equation (|55p becomes 



(a*+p) 



cPp dp f b 1 dp 



dfi^ d/i \ 6 2 rf/i 
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(59) 



To solve ((59)) . if we choose the initial epoch (^ = 
fj,o) to be a radiation dominated one (which is quite 
realistic in view of our current understanding of the 
universe) with c^ ~ 1/3, then from (|55|) we can 
easily see that c^ monotonically increases with /i. 
And in the interval (/zq < fJ. < oo) the function p{fi) 
is concave upwards. Therefore there must exist an 
earlier epoch at which p{p) > ji, which violates (0). 

4. Closed/open universe with varying sound speed: 
This is the most general case and let us try to 
find a solution with similar initial conditions as 
in the previous case. Since we know that very 
early universe was radiation dominated, let us 
suppose that there exists an epoch (ao << 1) with 



density /zq and pressure po where (hojPo) ~ l/^to- 
As we have already seen, ^R « l/og, hence the 
last term on the LHS of (|55l) becomes suppressed 
and in this case also one can easily show that c^ 
monotonically increases with /i. Therefore there 
must exist an earlier epoch ai < ao with pi > /io, 
where p{fJ.) > fi and ^ is violated. In other words, 
no solution satisfying © exists for (1551) that gives 
a radiation dominated era in the early universe. 

Hence for any physically realistic barotropic equation 
of state, if the new constraints are to be consistently 
propagated, we must have w°0 = 0. We thus proved an 
important theorem for shear- free perturbations of FLRW 
spacetimes: 

For an "almost" homogeneous and isotropic Universe 
filled with a barotropic perfect fluid subject to a physically 
realistic equation of state, if the fluid congruence is shear- 
free in a domain U, then it must be either vorticity-free 
or expansion-free in U. 

IV. DISCUSSIONS 

This result gives an interesting scenario. The linearised 
shearfree solutions - almost universally used to study the 
formation of structure by gravitational instability in the 
expanding universe, and believed to result in standard 
local Newtonian theory - do not have the same behaviour 
as shearfree Newtonian solutions. This may affect simple 
structure formation scenarios for rotating matter. 

Another interesting point that emerged from our anal- 
ysis is that there exists a class of barotropic equation of 
state (however unphysical that may be) for which the 
usual shear-free result can be avoided in the linearised 
case. It would be an interesting problem to see whether 
this same class of equations of state (or some similar 
class) allows shear free rotating and expanding solutions 
for the full non-linear Einstein equations for a barotropic 
perfect fluid. 
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